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We derive the pair-production probability in a constant electric field in Rindler 

O 

■ coordinates in a quasi-classical approximation. Our result is different from the pair- 

production probability in an inertial frame (Schwinger formula). In particular, it 
exhibits non-trivial dependence on rapidity and deviation from Gaussian behavior 

> : 

at small transverse momenta. Our results can be important for analysis of particle 



C3 



production in heavy-ion collisions. 



I. INTRODUCTION 



oo 
o 

o 

a,: 

Qh! In recent years a novel approach to particle production in high energy QCD has been 

developed. It is based on the observation that the color field of a high energy nucleus is a 



quasi-classical non-Abelian Weizsaecker- Williams field which strength increases with energy 

X ■ fi n 

(see reviews [l|, |2[ and references therein). The high field strength provides a semi-hard scale 
at which the coupling constant is small. This makes the problem of particle production in 
heavy ion collisions solvable in principle using the small coupling expansion with subsequent 
resumation of powers in strong field. This problem was solved at the lowest order in the 



density of color sources in 



4|. The resulting field is boost invariant, i.e. it depends on the 



invariant time r = v / 2x + x_ and in general, on transverse coordinate xj_, but not on pseudo- 



rapidity 7] = | In Moreover, at ear 
vanish while the longitudinal survive 



m 



y times after the collision r — > the transverse fields 
5j. Possible implications of this result were discussed 



6]. 



Existence of the strong longitudinal chromoelectric field in high energy nuclei collisions 
was first pointed out in the framework of the High Parton Density QCD in Q. In Ref. 8] 
it was further argued that dependence of the produced field potential on the transverse 



2 



coordinates can be neglected in the leading (transverse) logarithmic approximation implying 
vanishing longitudinal chromomagnetic field. 

The role of the strong longitudinal chromoelectricfield in particle production in heavy-ion 



collisions was discussed in many publications [20|, [21_|, [22|, [23|, [24], [27j . It was realized that 
the strong confining force produces particle pairs out of vacuum by means of the Schwinger 
mechanism [13I . . Since at high energies and in central collisions the nuclear color field 
is the Weizsaecker- Williams one, it became important to understand the pair production 
mechanism in the background of such field. The :>air production probability was calculated 
in the quasi-classical approximation in Refs. [7, 8|. It was argued in p| that the pair 
production is a mechanism by which the Color Glass Condensate can "evaporate" into 
thermalized quark-gluon medium over short time hence shading some light on a problem of 
fast thermalization in heavy-ion collisions at RHIC. 

In the present paper we reformulate the problem of pair production in heavy ion collisions 
in the Rindler coordinates which are the most natural for expressing the boost-invariance 



of the produced field 



AH, 
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12(| . In hydrodynamical models these coordinates are the 



coordinates of a fluid in a comoving frame. Since the Rindler coordinates correspond to the 
rest frame of the uniformly accelerated observer, the physical processes in the Rindler frame 
will in general look differently than in the inertial frame. In particular, we will argue that the 



probability of pair production by a constant chromoelectric field (IzTj) is significantly different 
from the Schwinger result [3] in contradiction to the common belief. This result may appear 
counterintuitive at first sight. However, I would like to stress that we not just use a different 
set of coordinates, but most importantly, quantize the system in a different space-time. It 
must be kept in mind that the Fock space in the Rindler space-time is physically different 
from the Fock space in the Minkowski space-time. In particular, the vacuum states \0)r and 
\0)m i n the Rindler and Minkowski space-time are different. This observation has been used 
by Unruh 25j to derive his famous effect: a uniformly accelerated detector (with acceleration 



a) becomes excited as if it was placed in a thermal 
other examples of that kind are considered in Ref. [ 



jath at temperature T = a/(27r). Many 
. We discuss this further in Sec. IVl 



The actual calculation of the pair production probability is carried out in the quasi- 



classical approximation pioneered in 



14 



16], see 17j for a review. 
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II. PAIR PRODUCTION IN THE RINDLER SPACE 

Consider a charged scalar field coupled to the Abelian gauge field A^. Equations of 
motion read 

D^D^cf) + m 2 <p = 0, (la) 
d^=f. (lb) 
where = <9 M + ieA^ and the current is 

j, = -i^D^ - (D^Yfl . (2) 

In absence of the back-reaction of a produced matter on the field the light-cone potentials 

n 

at the leading order in density of color charge read |3J 

A + = x + a(r) , A_ = — X-a{r) . (3) 

where we introduced the light-cone variables x± = (t ± x)/ \^2; t = y/2x + X- is the invariant 
time. 

Writing down we assumed that the field does not depend on transverse coordinates. 
Indeed, it was argued in [8] that production of low transverse momentum particles (pr Q s ) 
is dominated by the longitudinal chromo-electric fields. This implies that particles produced 
by longitudinal chromo-electic fields via the pair production mechanism constitute the bulk 
of the particle yield 2o| . 

In this approximation, the electric field is given by 

E z = --d T (ar 2 ). (4) 
r 

In this section we proceed neglecting the back-reaction and return to this problem later on 
in Sec. HVJ 

The Rindler coordinates (rj, r, xj_) are simply related to the Minkowski coordinates 
(t, x,x_i_) by the following relationships: 

r 2 = t 2 -* 2 , r^ln^. (5) 
2 £_ 

In the Rindler coordinates Eq. (flaj) reads 

-d T (Td T <p) - \dl$ - Vi0 + 2iead v (p + e 2 a 2 t 2 <p = 0. (6) 
r t 2 1 
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Switching over to the momentum space using 



yields the following equation for tp: 

(p + w(r)<p = 0, uj 2 (t) = -L + p l + fo. _ ear^j . (8) 

We can express the momentum p v in terms of the light-cone momenta of an inertial observer 

using the tensor transformation rule 

dx + dx_ , . 

Pr ' = ~~drf P ~ + ~drf P+ = X+P ~ ~ x ~ p+ - ^ 

For massless particles p v = 0. 

We are seeking solution of (jHJ) in the form = e lS . Introducing the auxiliary function 

u = S we have 

u 2 - zii - uj 2 = 0. (10) 

Assuming that S is slowly varying function so that \S\ < 5 2 we can adiabatically expand 
u = uo + «i + . . . where the functions uo and u\ satisfy — uj 2 = and 2uqU\ — iuo = 0. 
Thus, u = ± uj + iu/2uj + . . . and 

p_L,r) = _J_ e ±*/^, (11) 
'2u[t) 



which is a solution to (jHJ) in the WKB approximation. Using (jHJ) we can verify that the 
WKB approximation |^|/cu 2 <C 1 is valid at all r except at r — > 0. We will see below that 
the region r — > does not contribute to the pair production. 

In the constant field E = E z z we have according to (jlj) a — —E z /2. Assume for defi- 
niteness that a < 0, i. e. the field points out to the positive ^-direction. We have for the 
action 

S = J u(r)dr = J dr^l/Ar 2 + p\ + (pjr + e|a|r) 2 . (12) 

Changing the integration variable r' = e\a\r/p± and introducing a parameter a = e\a\/p 2 ± 
we can cast the integral in Eq. (fT2j) into form 



S 

where 



I J ^(t' 2 + A 2 + )(t» + A 2 _), (13) 



A l = \ ( 2 ^° + 1 ± y/(2p v a + l) 2 - a 2 (l + 4p^ 



(14) 
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Consider now two cases depending on the sign of the determinant D = (2p rj a + l) 2 — 
a 2 (l + 4p 2 ). Case a): D > 0. In this case the integrand of (fT3j) has a pole at r' = and 
two cuts along the imaginary axes — A + < Imr' < — v4_ < and < v4_ < Imr' < A + . In 
the WKB approximation the integration contour can be closed in the upper half plane [3] 
in which case the imaginary part equals 

= ~ [2p v a + 1 - ajl+4tf) (15b) 

where the Euclidean "time" is te = —it' . 

Case b) D < 0. A± become complex. It is convenient to change variables iz = r' 2 +p v a+^ 
so that the action reads 

S = -JL [ _? v /(b_ 2; )(b + 2; ) j 5 = i\fT)j2 G ReZ . (16) 

2a J iz - (p v a + |) 

The imaginary part of action arises from the discontinuity across the cut —B < Kez < B 
while the pole z = —i(p v a + 1/2) does not contribute. We have 

1 f° ^(B - z)(B + z) , 



71 

~Aa~ 



-2 Pri a-l + a^l+Ap2\ (17c) 
Combining the two cases (115b|) and (II Tel) we derive the final result 



lmS+ = - ^ | 2p„a + 1 - a yi+l^ | (18) 

Eq. (Tl8|) is not invariant under the transformation p,, — > — p n which implies by means of 
((7j) that cf) is not symmetric with respect to pseudo-rapidity inversion. This is an artifact 
of our choice of the electric field direction: a < means that E is pointing out in the 
direction of the positive z-axes (hence the subscript "+" of S in fflBl ). However, there is no 
preferred direction of the field since color charges are distributed approximately randomly 
in a nucleus. When the field points out to the negative z-direction the imaginary part of 
action reads: 



hnS- = - iL| - 2p v a + 1 - ayi + 4p 2 | (19) 
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The total probability of pair production is therefore 

w = ± (e~ 2Im5 + + e" 2ImS -) . (20) 

Note, that 2p v a+l — ayl + Ap 2 > when p v > (a 2 — l)/4a and negative otherwise, whereas 
—2p v a + 1 — a-^/I + 4p 2 > o when < (1 — a 2 ) /4a and negative otherwise. Therefore, we 
derive 

— ttItvJ , / 71 

w 



e _7rW cosh ^— (1 - ayjl + 4p 2 v )J , if > |a 2 - l|/4a . (21a) 

w = exp ~ a \j^- + %> 2 |) cosh(7r^) , if (pj < |a 2 - l|/4a . (21b) 

The pair production probability as a function of p v is demonstrated in Fig. [TJ 




FIG. 1: Rate of pair production w at a = 2 as a function of p^. 

It is clearly seen from Fig. [1] that the derivative of w with respect to p^ is divergent at 
the cusps located at p n = ±|a 2 — 1 1 /4a. Therefore, the quasi-classical approximation holds 
only in the three regions away from the cusps. 

For massless particles p v = and we have using (|21bl) 

w = , ^ > e£^ , (22a) 

w = e -f , < e E z . (22b) 
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The large p± tail of the probability distribution w, Eq. (122bl) . coincides with the Schwinger 
formula. However, in general, the spectrum of produced pairs in comoving coordinates is 
significantly different from the Schwinger formula even for massless particles (with p v — 0), 
see Fig. El 

1.2 
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FIG. 2: Rate of pair production w at p v = as a function of pt- The field strength is set to be 
E = 6/eGeV 2 so that Q 2 . ~ eE/ir ~ 2GeV 2 . Line marked by 'R' is the result is the Rindler-space 
calculation as given by whereas the line marked by 'S' is the Schwinger formula. Dashed part 
of 'R' corresponds to the values of pt at which the quasiclassical calculation breaks down. 

III. NON-ABELIAN CASE 

To generalize the obtained results for the case of SU(2) gauge interaction we write down 
equations of motion for the gauge field in the external field 

(-D 2 5,„ + 2ieF^)W, = 0, (23) 

where the gauge condition D^N 1 = has been used. Here we made a separation of the 
SU{2) gauge field into two parts: a quantized field and the classical background A^. 
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Eqs. ( 1231) can be diagonalized with respect to the light-cone gauge field components 
W± = W ± iW z 

- (D 2 ± 2eE g )W± = 0. (24) 

Solving these equations in the Rindler coordinates yields for W± the same equation as for 
0, see (IT]),® but with the shifted frequency uj 2 — > uj\ = u> 2 ± 2eE z . Therefore, the pair- 
production probability in the case of SU(2) can be read off the Eqs. ( 121 aft and ( I21bl) with a 
replaced by a±, where 

(25) 



eE K 



a± 



2(p 2 ± ±2eE z ) ' 

and sum over polarizations is performed in the arguments of cosh and exp: 



w 



e~ 2 ^ I cosh (V — I 



;i-a ±A /l + 4p, 



if > |4 - l|/4a± 



(26a) 



w = exp 



7T 



2a± 



(l_ a± /l + 4p 2) 



cosh(27rp^) , if < \a 2 ± — l\/4a± . (26b) 



We see that in non-Abelian case < a± < 1/4. As the result, asymptotic behavior of the 



pair-production probabilities is slightly modified with respect to the Abelian case (I22a[) . (I22bj) 
(for p v = 0) 

p 2 ± > eE z , (27a) 
pl«e^. (27b) 



w 



6 ei? 2 



W 



Large values of ImS* insure applicability of the quasi-classical approximation in these asymp- 
totic regions. 



IV. EFFECT OF BACK-REACTION 

To estimate the back-raection of produced particles on the background field note that the 
only non-vanishing light-cone components of the field-strength tensor are _F H = —F h = 
Eg, where E z = — i(? r (a(r)r 2 ). Therefore, Eqs. (Ilbj) can be written as 

dE z . 8E Z 
£ — j+ ; £ — j- _ (28) 

dx- ' dx + 

Switching to the Rindler coordinates r and r\ where pseudo-rapidity r\ = g hi — we find 
from flU that 

J+ = x +j( T ), J- =-X-j(r), (29) 
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which satisfies the current conservation = for arbitrary j(r), provided that j M is 

independent of the transverse coordinates. Finally, we obtain 

d T E z = -j{r)r. (30) 

The back-reaction effect can be neglected if \d T E z \ <C \E z \/t which implies that 

W)\ « ~ 2 \E Z \. (31) 

This is equivalent to requiring that the energy of interaction of the current with the external 
field j^A^ be much smaller than the energy of the field El/ 2 itself. 

To check the validity of the adiabatic approximation (13T1) we need to calculate the current 
induced by the external field. Substituting (CO) and ( ITTi) into (T5]) we derive in a classical limit 

(j(r)) = ^9 J ^fu{r) Or ~ 6aT ) 2w( ^' P±) ' (32) 

where g is the number of degrees of freedom of each of the produced particles. Formally, 
integral over p v diverges (since w is constant at large p v ). However, the rhs of fl32l) with 
w = 1 is just a vacuum term. It vanishes by reflection symmetry for the fixed integration 



3- 



boundaries of the kinetic momentum p v — ear 2 

In the physical situation we are interesting in (gluon production) p v w 0. Thus, we use 
Eqs. (127|) to estimate the induced current (|32l) . At later times r ^> pj 1 , (eE z )~ l integration 
over \p n \ < - and p\ G [0, oo) in fl32|) yields 

1 2eE z g * 

(j(r)) "^W 6 ■ (33) 
Then, employing (13"T1) we conclude that the back-reaction effect is small if the following 
condition is satisfied 

e g < 47r 3 e^ . (34) 
In the non-Abelian case the similar calculation results in even weaker condition 

e g < 87r 3 e 4 ^ , (35) 

which is satisfied for all reasonable values of coupling e and number of degrees of freedom 
g = 2{N 2 — 1). Note, that the average square of the transverse momentum of the produced 
particles is given by (p 2 ± ) = eE z /2Ti, see ( 127bl) . Therefore, the estimate (1351) is valid at 
r ^> l/^/2ir(p±), i.e. at all times of interest in heavy-ion collisions. 
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V. DISCUSSION 

The central result of our paper is Eqs. (I2T]) and (|26|) . They solve the problem of particle 
production by a constant electric field in the Rindler coordinates in absence of back-reaction 
of produced pairs on the field. We observed that the pair production probability in Rindler 
coordinates is essentially different from that observed by an inertial observer as given by the 
Schwinger formula, Fig. [2j Although the classical background field is boost invariant, the 
quantum corrections strongly violate it. 

It is important to emphasize, in order to avoid a possible confusion, that had we chosen 
to quantize the dynamical fields and in Minkowski coordinates and use the quasi- 
classical approach discussed in this paper, we would have obtained Schwinger's formula (see 
e. g. In the present paper, the dynamical fields are assumed to be quantized in the 

Rindler coordinates r and 77. At any given Rindler time r, the quasi-classical trajectories 
are = const, corresponding to the uniform acceleration in Minkowski space. On the 

contrary, Minkowski observers move along the straight lines. As a result, an observer in a 
Rindler frame is not equivalent to an observer in a Minkowski frame leading, in particular, 



to the famous Unruh effect 25j. In the present paper we argued that the pair production 
probability in a constant background field in Rindler space is also different from that in 
Minkowski space. 

The pair production probability w is proportional to the particle phase space density and 
can be used as an initial condition for the subsequent hydrodynamical evolution. It has 



been argued in [26| that small deviations from the boost invariance are amplified by the 
Weibel instability and can speed up the process of thermalization in high energy heavy-ion 
collisions. Therefore, our results can be important for understanding the final spectra of 
particles produced in heavy-ion collisions. 

Our estimate of back-reaction shows that the induced current is small. However, the 
quasi-classical approach to the back reaction is admittedly approximate. It ignores both 



the produced pairs correlations which may induce the background field oscillations 27| and 



quantum fluctuations which has to be properly renormalized 30j. We are going to address 
these important issues in forthcoming publications. 

It is also important to generalize the obtained results to include fermions. This will allow 
us to address the question of heavy quark production at small p± and shed light on the 
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heavy quark multiplicities. 
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